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Experimental  Observation  of  Fluid  Echoes  in  a 
Non-Neutral  Plasma 

Jonathan  H.  Yu  and  C.  Fred  Driscoll 


Dept,  of  Physics  and  Institute  for  Pure  and  Applied  Physical  Sciences, 

University  of  California  at  San  Diego,  La  Jolla  CA  92093-0319  USA 

Abstract.  Experimental  observation  of  a nonlinear  fluid  echo  is  presented  which  demonstrates 
the  reversible  nature  of  spatial  Landau  damping,  and  that  non-neutral  plasmas  behave  as  nearly 
ideal  2D  fluids.  These  experiments  are  performed  on  UCSD’s  CamV  Penning-Malmberg  trap  with 
magnetized  electron  plasmas.  An  initial  m,  = 2 diocotron  wave  is  excited,  and  the  received  wall 
signal  damps  away  in  about  5 wave  periods.  The  density  perturbation  filaments  are  observed  to 
wrap  up  as  the  wave  is  spatially  Landau  damped.  An  m,  = 4 “tickler”  wave  is  then  excited,  and  this 
wave  also  Landau  damps.  The  echo  consists  of  a spontaneous  appearance  of  a third  me  =2  wave 
after  the  responses  to  the  first  two  waves  have  inviscidly  damped  away.  The  appearance  time  of  the 
echo  agrees  with  theory,  and  data  suggests  the  echo  is  destroyed  at  least  partly  due  to  saturation. 


Longitudinal  plasma  wave  echoes  have  been  previously  observed  in  a neutral  plasma 
[1],  and  here  we  present  observation  of  diocotron  wave  echoes  in  a pure  electron  plasma. 
Diocotron  waves  vary  as  8 n(r,  8)  = 8n(r)exp[i(m8  - tot)],  with  essentially  no  z-variation. 
These  waves  can  be  considered  low  frequency  E x B drift  modes  representing  surface 
waves  or  bulk  shape  distortions  of  the  electron  column.  The  diocotron  mode  damping 
seen  here  is  a result  of  spatial  Landau  damping  [2].  The  echo’s  existence  depends  on  the 
reversibility  of  the  phase  mixing  process,  and  thus  demonstrates  the  reversible  nature  of 
spatial  Landau  damping. 

Theory  predicts  that  echoes  occur  in  ideal  (incompressible  and  inviscid)  2D  fluids  [3], 
and  the  echo  has  been  studied  numerically  in  2D  fluid  shear  flows  [4].  In  the  2D  E x B 
drift  approximation,  the  dynamics  of  electrons  in  a pure  electron  plasma  are  isomorphic 
to  that  of  an  ideal  2D  fluid.  The  density,  drift  velocity,  and  electrostatic  potential  of 
the  electron  plasma  correspond  to  the  vorticity,  velocity,  and  streamfunction  of  the  2D 
ideal  fluid,  respectively  [5].  Detecting  the  echo  in  a non-neutral  plasma  is  striking,  since 
one  might  think  that  non-ideal  effects  would  destroy  phase  information  in  the  density 
distribution. 

The  diocotron  waves  and  echoes  are  excited  and  received  using  wall  sectors  of  a 
cylindrical  Penning-Malmberg  trap.  Fig.  1 is  an  oscilloscope  trace  of  the  signal  received 
from  opposite  60°  wall  sectors.  This  detection  scheme  gives  a maximum  response  to 
m = 2 perturbations. 

An  initial  m,-  = 2 wave  is  excited  by  applying  V}  = 5 — 30V  to  two  opposite  60°  wall 
sectors  (different  from  the  detection  sectors)  for  about  one  half  of  an  m = 2 period.  The 
first  pulse  of  Fig.  1 is  the  signal  received  as  the  initial  wave  of  amplitude  A,-  decays  in 
about  five  oscillations  due  to  spatial  Landau  damping.  A second  mt  = 4 “tickler”  wave 
is  then  excited  a time  T after  the  initial  wave  is  excited.  The  tickler  wave  is  excited  by 
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FIGURE  1.  The  signal  received  from  two  opposite  sectors  of  a wall  cylinder  shows  the  initial,  tickler, 
and  echo  waves. 


applying  a square  pulse  of  V,  = 1 — 40V,  for  about  one  half  of  an  mt  = 4 period,  to 
four  30°  sectors  each  separated  by  90°.  The  second  pulse  of  Fig.  1 shows  the  time  at 
which  the  m,  =4  tickler  wave  is  excited.  This  signal  is  not  actually  the  m = 4 plasma 
oscillation,  but  rather  a direct  coupling  of  the  tickler  excitation  to  the  wall  detector.  The 
third  wave  packet  of  Fig.  1 is  the  received  echo,  and  the  peak  echo  response  occurs  at  a 
time  te  « 2x.  The  peak  echo  amplitude  is  denoted  Ae. 

We  have  also  imaged  the  spatial  evolution  of  the  electron  density  during  wave  damp- 
ing and  echo  generation.  At  a chosen  time  in  the  evolution,  the  plasma  is  dumped  onto  a 
phosphor  screen  and  recorded  by  a CCD  camera.  The  perturbation  images  shown  in  Fig. 
2 are  created  by  subtracting  a 0-symmetric  density  image  from  the  density  image  of  the 
perturbed  plasma,  giving  8n(r,0,r).  These  images  directly  represent  phase  space  evolu- 
tion, since  angular  momentum  is  proportional  to  r2  in  our  magnetized  electron  plasma. 
Thus,  the  filamentary  images  illustrate  the  theory  concepts  of  phase  mixing  (and  un- 
mixing). 

The  first  three  images  of  Fig.  2 show  the  initial  m,  = 2 wave  forming  spiral  filaments 
as  the  perturbation  evolves  in  the  sheared  background.  The  images  in  Fig.  2 correspond 
to  the  times  labelled  by  the  letters  in  Fig.  1 . The  m,  = 4 wave  is  excited  in  Fig.  2d, 
and  the  thin  filaments  are  the  remnants  of  the  initial  wave.  The  tickler  wave  nonlinearly 
interacts  with  these  remnants,  producing  an  me  = 2 echo  response  which  peaks  in  Fig. 
2h. 

Theory  gives  an  expression  for  the  time  of  the  echo  appearance,  and  this  expression 
agrees  with  experiments.  In  the  resonance  layer  where  0)£(r)  = (o/m , the  initial  pertur- 
bation can  be  written  as  8«,(r)exp[i>«,(0  - C0£-(r)/)].  Likewise,  the  tickler  perturbation 
varies  as  8n,(r)exp[(Vn,(0-  0)£(r)(t  - x))],  but  also  modulates  the  initial  wave.  This  cre- 
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FIGURE  2.  Density  perturbation  images  at  the  eight  successive  times  marked  in  Fig.  1 . Figs,  a,  d, 
and  h show  the  initial,  tickler,  and  echo  waves  at  maximal  amplitude.  Figs,  a - c show  the  spiral  density 
filamentation  and  wind-up  of  strong  spatial  Landau  damping.  Figs,  f - h show  the  “unwinding”  which 
creates  the  echo. 


ates  a second  order  perturbation  varying  as 

8ne(r,6,t)  — 8«,(r)8n,(r)exp[i(m,  - mj)Q+ im,0)£(r)x  - i(m,  — m,-)cO£(r)f].  (1) 

The  time  of  the  peak  echo  response  occurs  when  a radial  integral  over  the  perturbation 
in  Equation  (1)  does  not  phase  mix  to  zero,  and  this  occurs  when  the  phase  has  no  radial 
dependence  [6].  Equating  the  C0£  terms  of  Equation  (1)  to  zero  yields  the  time  of  the 
peak  echo  response 


te  = Xmt/{m,-mi).  (2) 

At  this  time  the  echo  perturbation  in  Equation  (1)  has  a 0 dependence  of  me  = mt  — rrij. 
This  mode  number  expression  agrees  with  the  observed  me  = 2 structure  of  the  echo. 

We  find  that  the  predicted  time  for  the  echo  response,  given  by  Equation  (2),  fits  the 
data.  Fig.  3 shows  the  measured  time  of  the  echo  response  te,  versus  the  time  of  the 
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FIGURE  3.  Measured  time  t , of  the  peak  echo  response  versus  tickler  delay  time  t. 


tickler  launch  x,  both  scaled  by  m — 2 wave  periods  72.  The  theoretical  line  te  = 2x  is 
shown  as  a dashed  line. 

The  initial  and  tickler  waves  do  not  have  symmetric  roles  in  their  interaction  leading 
to  the  echo.  Fig.  4 shows  the  echo  amplitude  Ae,  measured  at  t — te,  versus  x,  for  four 
different  (Vit  V,)  combinations.  Ae  increases  approximately  linearly  for  small  x,  and  then 
peaks  and  decreases  at  larger  X.  x^  is  defined  as  the  tickler  launch  time  which  results 
in  the  maximum  measured  echo  amplitude.  We  find  x^  does  not  depend  on  A,  (see 
dashed  curves  of  Fig.  4).  However,  the  solid  curves  of  Fig.  4 show  that  xmax  does  depend 
on  A,.  In  addition,  the  three  curves  with  V/  = 15V  give  nearly  the  same  Ae  at  Xmax- 

The  observed  asymmetry  in  initial  and  tickler  wave  roles  is  consistent  with  the  theory 
idea  of  saturation  destroying  the  echo  at  large  x.  A ballistic  theory  for  the  echo  amplitude 
yields  Ae(t  = te)  « A, J\  (pA,x),  where  J\  is  a Bessel  function  of  the  first  kind,  and  P is 
a function  of  m,  and  m,  [6],  A,  is  the  maximum  tickler  wave  amplitude  and  is  assumed 
to  be  proportional  to  the  applied  V,.  The  initial  wave  amplitude  A,  determines  the  echo 
amplitude  Ae  at  x^,  while  the  tickler  wave  amplitude  A,  determines  x^  itself.  These 
theory  results  agree  with  the  curves  in  Fig.  4,  indicating  that  saturation  is  at  least  partly 
responsible  for  the  decrease  of  Ae  at  large  x. 

Theory  of  echo  saturation  yields  as  many  as  300  m = 2 wave  periods  in  Xmax,  while 
experiments  give  no  more  than  about  50  m = 2 wave  periods.  The  discrepancy  between 
theory  and  experiment  at  small  tickler  wave  amplitudes  suggests  that  non-ideal  effects 
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FIGURE  4.  Measured  amplitude  Ae  of  the  echo  at  t = te  versus  T.  The  numbers  in  parenthesis  are  the 
initial  and  tickler  external  excitations  in  volts.  The  decrease  in  echo  amplitude  at  large  t is  due  at  least 
partly  to  saturation. 


may  destroy  the  echo.  One  non-ideal  effect  is  the  magnetron  drifts  resulting  from  the 
end  confinement  fields,  which  cause  9-smearing  [7]  and  which  would  presumably  kill 
the  echo.  In  addition  to  9-smearing,  viscosity  can  also  wreak  havoc  on  the  ideal  2D  fluid 
isomorphism.  As  the  density  perturbation  filaments  wind  up,  the  radial  spacing  between 
filaments  decreases  and  the  filaments  become  thinner.  After  sufficient  spiral  wind-up 
the  radial  gradients  would  approach  a magnitude  for  which  viscosity  is  appreciable,  and 
phase  information  would  be  destroyed.  The  influence  of  these  non-ideal  effects  on  the 
echo  will  be  pursued  in  future  work. 
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